ADA100740 


lEVB^  d 


T) 


^  JiTATXSTICAL  PROPERTIES  OF  A  JRANDOM  ARRAY 
OF  ACOUSTIC  SENSORS  IN  A  MULTIPATH  ^ENVIRONMEN^ 

i  ii)  )  Fred /Haber  \ 


^  juHaa'9* 


r//i 


9  ‘  /’ 


’■-.in.'jal  Libur 


JPraov.-:i)  PJTOxi  Mr-Eiig* 


^  'ij 

Valley  Forge  Research  Center 
The  Moore  School  of  Electrical  Engineering 
Uni.versity  of  Pennsylvania 
Philadeli/hia,  Pentisiylvania 


'  i9b'U 


II  I.  .  .^v 


hV-  f  ■/ 1‘‘ 


;  /-■,  • 


ABSTRACT 


r/J  ^  - 

iim  and  main  beam  gain 


The  statistical  properties  of  the  power  patteim  and  main  bean  gain 
are  here  determined  for  an  array  of  randomly  located  submerged  acoustic 
sensors.  Tl.e  system  Investigated  models  a  sparse  array  in  which  the  sen- 
Gors  are  dropped  haphazardly  over  a  region  or  in  which  the  sensors,  how¬ 
ever  they  arc  initially  placed,  become  spatially  diffused  by  a  process 
akin  to  a  two-dimensional  random  walk.  Signal  energy  is  assumed  to  arrive 
over  a  vertically  dispersive  channel  typical  of  the  long  range  deep  sea 
acoustic  channel  with  a  blgradient  sound  speed  profile.  Results  are  ob¬ 
tained  for  the  mean  value  of  the  power  pattern  and  the  mean  and  variance 
of  the  main  beam  pov.cr  gain  as  a  function  of  arrsy  size.  It  is  shown  that 
in  typical  cases  the  mean  array  gain  will  preserve  its  value  within  3  dB 
until  the  dlspersiou  parameter  of  the  array  measured  by  the  standard  dev¬ 
iation  of  element  location  is  around  35  wavelengths.  The  results  are  fouvid 
to  be  consistent  with  results  obt.uined  by  others  for  the  coherence  distance 
in  a  multipath  acoustic  field, 
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STAJISTICAL  PROPERTIES  OF  A  RANDOM  ARRAY  OF 
ACOUSTIC  SENSORS  IN  A  MULTIPATH  ENVIRONMENT 

INTRODUCTION 

The  statistical  properties  of  the  power  pattern  of  an  array  of  acous¬ 
tic  sensors  suspended  froa  individual  freely  floating  buoys «  receiving  froa 
a  source  via  a  tiae-disperslve  aedlua  are  here  investigated.  '  The  systea 
geoaetry  is  indicated  in  Figure  1.  The  array  elements »  numbering  N,  are 
assumed^  distributed  in  a  region  centered  on  the  origin  of  coordinates.  By 
independent  aeans  the  systea  learns  the  position  y^»  ,  1  -  1,2,  ...  N 

of  each  of  its  elements.  We  suppose  that  the  array  then  organizes  Itself 
at  a  specified  frequency  forming  a  beam  aimed  in  some  selected  direction 
by  adding  suitably  phased  versions  of  the  element  outputs.  If,  for  Instance, 
the  element  outputs  were  of  unit  magnitude  at  the  specified  frequency,  if 
the  signal  wert  arriving  at  angle  (9,i^),  and  if  the  array  were  focused  to 


■  ■  I,  t 
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receive  a  signal  Tay>  arriving  in  the  plane  at  a  colatitude  angle  6  the 

8 

conplex  array  output  at  that  frequency  would  be, 


N  jk[]t  slndcos^  4-  y  (slnOsln^ 

A-  X  e  “ 

n-1  ; 

1 


sine  )  +  z^Ccosft  -  C086  )1 
8  n  8 


(1) 


tffaere  the  wavenumber  k  ■  2Tr/X*  X  being  the  wavelength.  In 
this  work  we  will  assume  that  all  elements  are  at  the  same  depth,  all  in 
the  x-y  plane,  so  that  all  0.  Though  surface  waves  will  cause  verti¬ 

cal  displacement  of  the  elements,  the  wavelengths  of  interest  are  such  that 
« 

in  placid  seas  (sea  state  <  4)  the  displacement  is  less  than  0.1  wavelength. 
Furthermore,  the  system  envisioned  is  expect'^d  to  utilize  a  drouge  with 
each  buoyed  element  so  that  vertical  motion  will  be  filtered.  Element  po¬ 
sitions  in  the  (x,y)  plane  will  be  assumed  independently  distributed  accord¬ 
ing  to  some  appropriate  two-dlmentsional  probability  density  function  (pdf) 
£^(x,y),  as  will  be  discussed  later. 

The  dispersion  model  utilized  assumes  a  collection  of  M  planar  wave- 
fronts  impinging  on  the  array  all  originating  from  the  same  source  and 
all  arriving  with  the  same  azimuth  angle  but  X7lth  different  colaMtude  an¬ 
gles.  This  model  requires  that  the  phase  front  corresponding  to  a  given 
ray  arriving  at  the  array  center  be  adequately  approximated  by  a  planar 
surface  wherever  the  phase  front  contacts  the  array.  If  the  array  length 
in  the  azimuthal  direction  of  arrival  of  the  ray  is  d  and  the  colatitude 


angle  of  the  ray  is  6  the  distance 


across  the  wave  front  over  which  planarity 


should  hold  is  dcosO.  Each  wivefront  is  characterized  by  a  complex  amplitude 


B  e 
m 


ffl 


m**  It^  ...  M, 


(2) 
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at  the  frequency,  frthe  phases  and  amplitudes  are  measured  at  the  origin 


of  coordinates,  the  former  relative  to  an  arbitrary 
be  treated  as  a  random  variable.  Independent  of 
dom  variable  uniformly  distributed  In  2'n'  and  ()>  for 

IB 

Burned  Independent.  .No  assumption  is. made  about  the 
for  different  m.  The  colatitude  angle  of  arrival 


reference.  B  will 
m 

d)  will  be  assumed  a  ran- 
different  m  will  be  as- 
dependence  among  the  S 

B 

\  , 

of  a  wavefront  will  be 


denoted  6  ,  m*  1,2,  ...  M.  The  6  will  be  viewed  as  nonrandom  constants. 

®  o  .... 

They  may  be  taken  to  be  equally  spaced  angular  samples.  The  modelj employed 
corresponds  to  one  used  by  Smith  [Ij  to  calculate  spatial  coherence  In  a 
multipath  channel. 

With  the  wavefronts  arriving  at  azimuth  angle  ^  and  the  array  focused 
to  receive  a  plane  wave  from  a  source  at  azimuth  angle  —  ,  l.e.  from  a 
source  in  the  j—z  plane,  and  colatitude  angle  6^  the  total  array  output  bi- 
comes 


M  N  Jk[x  sln6  cosi^  +  y  (slnB  sin^  -  8ln9_)  +  ♦  ] 

A/x  n  m  8  m 

A(*,e3)  -  I  I  B^e 

This  Is  the  conplex  array  pattern.  The  statistical  properties  of  the  cor- 

I  2  I  * 

responding  array  power  pattern  |A  ($,63)1  will  be  investigated  below. 


The  terminology  used  here  differs  somewhat  from  that  given  In  Urlck  (2], 

We  Initially  calculate  an  array  pattern  |A^(i;’,0  )|  which  corresponds  to 
the  square  of  Uiick *3  response  function  [l,pp.  A9-50].  We  then  de¬ 

termine  a  normalized  mean  array  pattern  which  Is  slr.llar  to  Urlck's  "beam 
pattern"  [2,  p.  50].  The  normalization  used  here  differs  however  from  that 
used  by  Urlck  allowing  us  to  account  for  loss  of  coherence  across  the  array. 
Finally,  we  define  the  normalized  mean  array  pattern  evaluated  on  the  main 
beam  as  the  mean  power  gain. 

[1]  P.  W.  Smith,  Jr.,  "Spntlnl  Coherence  in  Multipath  or  Multimodal  Chan¬ 
nels,"  Jour.  Acoustic.  Society  of  America,  Vol.  60,  No,  2,  August  1976, 
pp.  305-310. 

[2]  R.  J.  UiXck,  Prlnctplen  of  UnJervater  Sound,  McOraw  Hill  Book  Co., 
Second  Edition,  1973. 


STATISTICAL  PROPERTIES  OF  ARRAY  POWER  PATTERN 


nie  mean  value  of  the  array  power  pattern  is  given  by 


M  M  M  M 


I  I  I<.  B  > 

nj^“l  ^  ^ 

jklx^  co8(^-x^  co8<|)+yj^^(ain0^^8in<J)“8ineg)~yj^^(sih0jj^ein4)“8ln0^)f;;>> 


expectation  <Cexp(ji^^  -  ^  when  nij^  »•  02  zero  other  ¬ 


wise  so  that  (A)  l8 


M  N  N  ,  jk(x  -X  )8in0  cos* 

<:lA^*,e.)l>-  I  I  I  <=^><*  1  ^ 

n-l  n,"l  n»«l 


‘1  *  “2 


Jk(yn  )  (slnOj^^sln*  -  slu'i^)^^ 

^  s  12 

(5) 

Using  the  assumption  that  the  element  positions  are  Independent  random  vec¬ 


tors  (5)  Is  written 


M  -  N  N  jk(x  -jx  )ein0  C08* 

<.|A^*.e,)  |>.  I  <B>  IH  .1  I  "f  B 

in“l  n,"l  no"i 


1-x  Mj- 
Hjllnj 


1  “  ’ 

B-1  In-1 


jk(y  -y  )  (slnG^sln*  -  8in0^)^> 
•  c  ”1  *^2 

j2 

,  N  jklx  sine  cos*  +  y  (slnO  sin*  -  slnO  )1^ 

+  I<“  1 


*  N  Jk[x  sine  cos*  +  y  (sinG  aln*  -  sinO  )]  _ 

-  I  n  m  ®i>r  (6) 
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The  expectations  Inside  the  brackets  are  two  dimensional  characteristic |func- 
tlc-ns  of  tht?  random  vectors  Assuming  all  elements  to  have  Identl- 

C'jilly  distributed  location  vectors  and  denoting 


^  Jk[x  sln0  cos^  +  y^(slo5i  slni^  -  8ln0^)_  I 

<e  "  "  “  “  I 


then 


M 


m«l  ' 


(8) 


At  this  point  we  specialize  the  distribution  of  the  location  vectors. 
We  assume  the  effect  of  the  forces  tending  to  scatter  the  array  elements 
to  be  modeled  by  a  two  dimensional  random  walk  with  Independent  increments 
along  the  coordinate  axes.  will,  after  a  time,  be  distributed  ac¬ 


cording  to  a  two  dimensional  random  variable  approaching  a  nonuil  wlthj  vari- 

2  2 

ance  along  x  and  y  given  by  o  and  a  respectively.  In  this  case 

X  y 


♦  (0  ,^,0  )  « 

xy'  m*^*  e\  , 


-[O^^sln^e  cos^4>  +  a^k^(sin0  siiufr  -  8lne^)^]/2 
g  X  n  y  a  ^ 


(9) 


The  model  chosen  accounts  for  element  diffusion  but  Ignores  drift  compon¬ 
ents  which  are  sure  to  be  present.  The  assumption  Is  Implied  that  trans¬ 
lation  of  the  entire  .<.rray  will  not  seriously  affect  array  response  when 

attempting  to  focus  on  a  distant  target  If  \  he  translation  Is  small.  Wc 

2'  2 

point  out  that  a  and  a  are  functions  of  time;  as  the  array  ages  these 
X  y 

parameters,  which  are  a  measure  of  the  size  of  the  array,  will  grow. 


-6- 


I£  we  were  to  suppose  that  initially  the  elements  are  close  to  one 
another  so  that 


«  (6  )  -  1 
xy  m*  ’  8 


)1>-  I  CBf>  -  nV  (11) 

■  J  ^  ^  Is  the  mean  square  value  of  the  total  signal  arriving  at  the 

array  center  and  (11)  represents  the  power  delivered  by  the  array  when  it 
is  sufficiently  small  not  to  be  defocused  by  the  multipath.'  It  is  useful 

to  normalise  lA^(<ji,e  )|  by  We  denote  this  random  varlable»' 

8 


r  -  -■ 


|A^<fr,e  )| 


tlie  normalized  array  power  gain.  Its  mean  value  is 


14.  n  -oVsin^e^cos^*  -aV(slne^8in^-slne^)- 

^  N^il  b2  *  (13) 


*An  alternative  normalization  is  given  by 

lA^(+,e^)l 

**1  "  -  M  j*  o 

-'I  I  »„«  ”1' 

®®*1 

Thl.  ra.0o„  variable 

to  small  size.  T  as  given  by  (12;  pp  total  arriving  signal 

1.  unity,  in  ‘he  variability  asa.Ua«^^^  „|  ,,, 

power  at  the  array  center  has  J ^  ^^ray  defocuslng. 

quantity  is  more  nearly  representative  ^  there- 

r.  la  howavur  a  morn  difficult  auantity  wltl.  which  « 
fire  settled  on  T  In  which  the  nornllimtlon  la  done  witlv  e  constant. 
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As  a  final  step  we  define 


B(e  )Ae 

n 


(14) 


that  Is,  we  suppose  the  continuum  of  possible  signal  arrival  angles  to  be 
quantized  Into  Increments  A6  and  that  the  fractional  power  obtained  from 
the  mth  Increment  Is  B(8^)A6.  By  allowing  the  Increments  to  become  small 
(13)  will  be  approximated  by  cr.  Integral  as  follows: 

1  1  -J  \^8ln^0cos^4>  “O^lt^(8ln0sln(^-sln6  )^ 

<r>  -  ~+ (1-i)  /  B(e)e  y  .  *  d0 

I  (15) 


<r>  as  given  by  (IS)  Is  shown  evaluated  In  Figure  2  for  the  case 

0 

2  2  2  I 

o*  ■  o*  ■  0  wtth  ok/2iT  “  a/X,  a  family  parameter,  given  by  5,10,20,  and  40  wave' 

^  #  I  , 

lengths.  The  power  density  3(6)  Is  assumed  uniformly  distributed  over  ilO* 
relative  to  the  horizontal.  The  array  Is  assumed  focused  for  a  source  in 
the  plane  of  the  array,  ^that  is  0^  ■■  90^,  and  at  an  azimuth  angle  of  90*^. 

Of  particular  Interest  Is  the  magnitude  of  the  normalized  mean  array 
pattern  evaluated  on  the  main  beam  as  a  function  of  array  size.  We  refer 
to  this  quantity  as  the  power  gain  Its  mean  value,  la  obtained 

from  (13)  or  (IS)  evaluated  at  ^  90*.  Thus  using  the  discrete  ray  model 


(16) 


FIGURE  2 


.  PATTERN  OF  NORMALIZED  MEAN 


ARRAY 

POWER  GAIN. 

(a) 

o/A  • 

5 

(b) 

o/A  ■ 

10 

(c) 

o/A  " 

20 

(d) 

o/A  - 

AO 

I 

Or,  using  the  continuous  approxlnaticn  for  the  integral , 

,  ,  -0^^(sln6-sln6^)^ 

<rr^-f  +  a~)  /  0(e>e  ®  dS  (17) 

In  (16)  and  (17)  the  subscript  on  has  been  dropped.  (17)  has  been 
evaluated  numerically  as  a  function  of  the  normalized  size  variable 

ok/2ff  «  c/X,  with  the  angular  distribution  of  energy  arriving^  6(6),  as  a 
parameter  and  6^  set  to  zero.\^  8(8)  was  assumed  uniform  over  angles 
is*,  ilO*,  and  i20X  relative  to  the  plane  of  the  array.  Results  are  shown 
In  Figure  3,  For  distant  sources  the  arrival  angles  are  apt  to  be  within 

NORMALIZED  MEAN  POWER  GAIN 


FIGURE  3. 1  MEAN  POWER  GAIN  AS  A  FUNCTION  OF  NORMALIZED 
j  ARTAY  SIZE  VARIABLE,  o/X. 

(a)  8(0)  1/10,  85*  <6<  95* 

(b)  8(6)  -  1/20,  80*  <B<  100’ 

(c)  8(0)  -  1/40,  70*  <e<  110* 


-,%r 


ilO'*  .  Note  that  delivered  power  Is  reduced  to  about  1/2  when  o/X  *  35  wave¬ 
lengths  J  For  sources  nearby,  bottom  and  top  reflections  may  result  In 
ene.-gy  arriving  at  steeper  angles  and  the  +20%  distribution  may  be  viewed 
as  a  model  suggesting  the  effect  In  such  a  case.  Here  delivered  power  is 
reduced  to  about  1/2  where  o/X  *  8  wavelengths. 

Setting  the  angle  6^  to  aero  means  focusing  the  array  for.  signal  ar¬ 
rivals  in  the  plane  of  the  array.  This  is  not  optimum  for  signals  arriv¬ 
ing  over  a  dispersion  of  latitude  angles.  To  show  this  we  have  plotted 

this  mean  power  gain  as  a  function  of  the  colatitude  angle  6  for  the  case 

s 

of  8(6)  uniform  in  +10*  around  the  horizontal  and  o/X  -  33.3  wavelengths. 

The  result  Is  shown  in  Figure  4.  Tlte  maximum  is  seen  to  occur  with  the 

beam  formed  for  6  ■  86*  (by  symmetry  it  will  also  bo  maximum  for  8  ■  94*). 

■  s 


80  8Z  84  66  88  90 

COLATITUDE  ANGLE,  0, 


FIGURE  4.  MEAN  POWER  GAIN  AS  A  FUNCTION  OF 
VERTICAL  AIMING  ANGLE  0^. 

B(e)  -  1/20,  80*  <8<  100“ 


:  -J 
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The  variance  o£>  the  power  gain  will  be  useful  as  an  Indicator  of  the 
gain  variability.  We  thus  determine  Var  )  where 

r  ■  |a^  (■?’,  0  )  1/N^B^.  Starting  f3)  evaluated  at  ^  ■  Tr/2  we  have 

O  M  8 


^  "  I,  X  '®ta' 

in«l  n*!! 


M  N  J[ky^(9in0^  -  sinO^)  +  <|>^1 


The  fourth  moment  of  the  magnitude  is  given  by 


C|A^1>>-  I  I  -<[B  B  b  B  ^ 

o'*^  ^  *  m,  m..  m, 

m^ym^tTDy  n^tn2in2»  12  3  4 


-  1  n^  -  1 


^i(k  y  -k  y  +k  y  -k  y 
.<i  “l  "l  "‘2  “2  ”3  "3  ”4  "4 


Me  have  used  the  abbreviated  notation  k  ■  kCsinO  -sinO^).  With  the  4 

m^  8  in^ 

Independent  uniformly  distributed  random  variables  in  (0,27r)|(19)  reduces  tc 


«  ,  „  Jlk  (y_  -y„  (y_  -7-  )] _ 

i  N  f  M  0%  ^  n  ®i  ^1  ®'a  ®‘!l  “A 

■  Z  I  Z  ^  ^  ^  3  3  4  > 


®l*V^  1  2 

“ij’*  ’«2 


J  [k  (y  -y  )  +  k  (y  -y  )  ] 

^3 


M  (y  -y  -^y  -y 

*  Z<»><e  '  ' 

m*l 
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SummaMon  over  the  n^,  1  »  l,2,3,4.,ls  now  carried  out  using  the  assump¬ 
tion  that  the  positions  of  the  different  elements,  the  )c-l,2,  ...N, 
are  independent.  The  summation  Is  straightforward  although  laborious. 
Carrying  out  the  nteps  we  can  then  write 


Var  r  » 
o 


<ia;i>  -  <|a;i>^ 

1I^b‘ 


H  M  <}l><<>  r  -“'“Li 

tBj“l  in2"l  B  ^  i. 

r-a^(k  •^k 

+  2N(N-1)  |e  ^+c  J 

2  m.  ■  m,  2  m,  m-  2  m.  in2 

+  4N(N-l)(N-2)c  ^^+e  J 


L2N(N-l)(N-2)(N-3)-  N^(N-l)^i  e 


-o^(k^  +k^  ^ 

“1  *“2  f 


+  I  - - - S -  .  -L  'jN(2N-1)  +  4N(N-l)^e 


2  ^  m 


+  N(N-l)(N-2)(N-3)e 


-2o^k^ 

m 


+  2N(N-l)(N-2)  e 


-3o^k^ 

m 


+  N(N-l)  e 


-4a^k^l 


The  ratio  ’  r  F  /'^V  where  the  numerator  is  given  by  (21)  and 

o  0 

the  denoiainatuA  ^,y  the  square  of  the  mean  power  gain  given  in  (16),  Is 
a  useful  measure  of  relative  variance.  The  result  Is  cumbersome. 


K  -  -  -"iuMf, 
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however,  and  expressions  applicable  to  j.ii"tting  cases  are  instructive. 

« 

Tv70  such  cases  are  here  evaluated  assuming  the  ray  amplitudes,  B  ,  are 

k 

equal  for  all  m  and  constant  so  that <11.2  .-/B  =>  (1/M) 2,  k  even. 

m  *■ 

The  first  case  treated  assumes  the  elemeuts  very  widely  dispersed 
2  2 

so  that  Ok  »  1  for  all  rays  except  ones  for  which  sin0  ■  sir9„. 

nx  ns 

Since  k  ••  k  (sln6  -  sinO  )  a  ray  along  the  aiming  angle  will  result  in 
m  A  8 

■  0.  Assuming  one  ray  is  along  the  aiming  angle  set  at  6^  ■  90° 
we  get,  using  (16)  and  (21), 


Var  r 
_ < 

<ro>' 


+  M(2W^-2M-1)  -  (2N^--N-1) 
2 


N(M+N-1) 


(22) 


Vfhen  M  “  1,  there  is  only  1  ray,  and  that  one  along  the  aiming  angle. 
There  is  no  multipath  and  the  array  will  be  correctly  focused.  The 
ratio  above  is  then  zero.  When  M  gets  large  without  bound  while  N  re¬ 
mains  finite  the  ratio  approaches  unity.  This  result  can  be  anticipated. 
For  M  large  the  elements,  being  widely  dispersed,  see  a  sinusoid  with 
Rayleigh  magnitude  and  random  phase.  The  ratio  in  (22)  in  that  case  is 
that  of  the  variance  and  squared  mean  of  an  exponential  random  variable 
for  which  this  ratio  is  unity.  If  M  is  held  finite  while  M  is  allowed 
to  increase  without  bound  the  ratio  tends  to  zero.  This  result  arises 
because  the  element  outputs  caused  by  the  one  ray  along  the  aiming 
angle  are  coherently  combined  by  the  array.  The  random  component  con¬ 
tributed  by  the  rays  off-axis  add  up  non-cbhcrcntly  at  the  array  output. 
The  latter  are  the  variance  producing  components.  But  as  M  Increases 


This  result  assumes  also  that  slnk^  /slnkn,  for  all  m./m,.  Should  there 

12  X  z 

be  rays  arriving  symmetrically  relative  to  the  horizontal  there  will  be 
some  slnk^  “sinkj^  for  m.^mT.  In  such  a  case  some  additional  terms  will 
be  required  from  ^21)  and  (22)  may  exceed  unity. 
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The  ratio  of  the  non-coherent  components  tends  to  zero. 

We  point  out  that,  if  no  ray  comes  in  at  an  angle  sufficiently  close 

2 

to  0^  to  make  a  significant  coherent  contribution  then  Vat 
approaches  (l-l/MN).  Now  the  ratio  approaches  unity  with  increasing 
M  or  N  as  one  would  expect. 


The  second  case  treated  is  one  for  which  the  element  locations  are 

reasonably  compact,  say  within  o/X  -  10,  and  N  is  large.  In  this  case, 

2  2 

if  ray  arrivals  are  within  +  10*  of  the  vertical  focusing  angle  o  k 
is  small  and  e  is  close  to  unity  (it  is  0.748  for  0  ■  10*).  The 

predominant  terms  in  (21)  are  then  those  of  highest  degree  in  N.  Extract- 


i 


In  (23)  we  have  set  such  factors  as  (N-1) (N-2) (N-3)/N^  to  unity  in  con¬ 
formity  with  specification  that  N  be  large.  Again  now  with  a  constant 
and  equal  for  all  n 
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2 

is  obtained  from  (16)  uitlng  the  assumptlone  pertinent  to  thle 
case.  The  ratio  Inside  the  bracketed  facto.  In  (24)  Is  1/M  for  0*0. 

For  a  ••  0  (24)  Is  Identically  the  variance  of  the  squared  amplitude  of 
the  sum  of  M  equal  amplitude  slnasoids  with  independent  phases,  all 
uniformly  distributed  In  27r;  this  reeult  for  a  0  simply  reflects 
the  variability  of  the  Incoming  total  signal  magnitude  at  the  array 
center.  It  is  Interesting  to  observe  that  If  B  were  Rayleigh  dls- 
trlbuted,  Var  would  be  unity  for  all  M  and  for  all  o/X 

for  which  (23)  Is  valid. 

To  convey  some  Idea  of  how  the  variance  changes  with  the  array  size 

parameter  a/X,  and  to  provide  some  results  without  approxltrntions  the 

2 

ratio  Var  was  calculated  for  some  representative  cases 

using  (16)  and  (21)  as  they  stand.  Table  1  shows  these  results  for  the 

case  of  M  rays,  M  ■  3,  5,  11,  and  21,  the  rays  arriving  at  equally 

spaced  angles  In  an  interval  of  t.  10°  relative  to  the  horizontal.  The 

beam  pointing  angle  was  set  along  the  horizontal  (6^  90^)  and  in  the 

direction  of  the  source.  The  number  of  sensors  N  was  taken  to  be  31. 

It  Is  worth  noting  that  for  M  •  3,  N  -  31  In  the  limiting  expression 

(22)  (o/X-H*),  Var  V  f<dV  "  0.1175  only  slightly  below  the  calculat* 

o  0 

ed  value  In  Table  1  for  M  "  3,  and  a/X  60. 

There  may  be  another  phenomenon  responsible  for  the  trend  of  the 
calculated  values  In  Table  1.  With  M  >  1  the  signal  power  at  the  array 
center  is  a  random  variable  as  a  consequence  of  the  random  phases  of  the 
incoming  rays.  The  ratio  of  Var  Is,  as  pointed  out  below 

(24),  given  by  (1-1/M).  This  is  exactly  the  value  in  Table  1  for 


>'\Q 


o/X  ■  0.  As  o/X  increases  the  ratio  Var  ^  decr’ases,  partly 
for  the  reason  discussed  under  (24) •  and  perhaps  also  because  the  signa 
amplitude  which  varies  with  position  is  averaged  by  the  array. 


M 


3 

5 

11 

21 

0 

0.67 

0.80 

0.91 

0.95 

5 

0.66 

0.80 

0.91 

0.95 

10 

0.62 

0.78 

0.91 

0.95 

15 

0.42 

0.73 

0.90 

0.95 

20 

0.20 

0.70 

0.88 

0.95 

40 

0.12 

0.64 

0.84 

0.92 

60 

0.12 

0.47 

0.81 

0.91 

TABLE  1.  RATIO  OF  VARIANCE  TO  SQUARED  MEAN  OF  ARRAY 
OUTPUT  AS  A  FUNCTION  OF  ARRAY  SI.’.E,  o/X,  AND  NUMBER, 
M,  OF  EQUAL  AMPLITUDE  RAYS,  RAYS  ARE  ASSUMED  EQUALLY 
SPACED  IN  AN  INTERVAL  ±  10  FROM  HORIZONTAL. 
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CONCLUSION 

Mean  po'ver  pattern,  moan  main  beam  gain  and  main  beam  gain  variance 
have  been  determined  for  an  acoustic  array  of  widely  scattered  submerged 
elements.  The  elements  are  assumed  to  be  organized  to  accept  a  plane 
wave  but  the  array  sees  a  multipath  field  typical  of  the  acoustic  field 
at  a  great  distance  from  a  source  In  the  deep  ocean.  The  results  obtained 
show  the  diminution  of  array  effectiveness  as  the  array  size  grows;  for 
a  typical  case  where  the  range  of  latitude  angles  of  the  arriving  signal 
Is  10*  above  and  below  the  horizontal,  the  mean  power  gain  fal!*  **.  3  dB  when 
the  element  spread  as  measured  by  the  element  position  standard  devia¬ 
tion  Is  35  wavelengths. 

In  earlier  work  Smith  [1]  calculated  the  normalized  coherence  magni¬ 
tude  as  a  function  of  horizontal  separation  in  a  long-range  transmission 
channel  with  bigradlent  sound  speed  profile.  In  particular  he  explicitly 
obtains  the  coherence  distance  for  50%  coherence  for  a  source  on  axis  and 
a  receiver  close  to  the  channel  edge.  It  turns  out  to  be  46  Wi.velengths 
If  the  difference  In  sound  speed  between  the  axis  and  the  receiver  loca¬ 
tion  is  20m/sec.  The  range  of  vertical  angles  of  arrival  is  dr  9.36*  rela¬ 
tive  to  the  horizontal  at  the  receiver  and  the  energy  density,  obtained 
from  earlier  results,  Is  assumed  by  him  to  be  uniform  over  the  range  of 
arrival  angles.  This  situation  Is  roughly  the  same  as  that  used  to  ar¬ 
rive  at  the  35  wavelength  figure  mentioned  above.  The  two  results  are 
interestingly  similar.  This  ought  to  come  as  no  surprise;  since  the  co¬ 
herence  distance  corresponds  to  the  array  dimension  useful  for  coherent 
comblnatloii  of  the  spatially  sampled  field.  The  crucial  factor  In  de¬ 
termining  the  coherence  distance  is  the  range  of  angles  of  arrival.  To 
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see  this  intuitively',  imagine  two  multipath  rays,  one  horizontal  and  the 
other  smaller  in  amplitude  and  at  6*  to  the  horizontal.  The  phase  ij*  of 
the  resultant  vs  horizontal  distance,  x,  is  expressible  as  [3] 

\(»(x)  -  kx  +  ♦(x) 

where  k  is  the  wave  number,  ^(x),  is  periodic  with  period  2‘ir/k(l*co86) 

X/ (1-0080).  With  6  ■  10*  the  period  is  about  66X.  This  is  not  too  far 
different  from  the  2o  range  (-  70  wavelengths)  for  an  array  power  gain  of 
O.S,  or  from  the  46  wavelengths  obtained  by  Smith  for  coherence  distance 
for  O.S  normalized  coherence. 


[3]  F.  Haber,  "Phase  Voriatlons  with  Position  in  an,  Underwater  Multi 
path  Fnvironinent  and  its  Effect  on  Array  Pattern,"  Valley  Forge 
Research  Center  Quarterly  Progress  Report  No.  2A,  Unl-verslty  of 
Pennsylvania,  the  Moore  School  of  Electrical  Engineering,  Phila¬ 
delphia.  PA  1910A,  pp.  20-30. 
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